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Abstract 

We show that active transport processes in biological systems can 
be understood through a local equilibrium description formulated at the 
mesoscale, the scale to describe stochastic processes. This new approach 
uses the method established by nonequilibrium thermodynamics to ac- 
count for the irreversible processes occurring at this scale and provides 
nonlinear kinetic equations for the rates in terms of the driving forces. 
The results show that the application domain of nonequilibrium thermo- 
dynamics method to biological systems goes beyond the linear domain. A 
model for transport of Ca^+ by the Ca^+-ATPase, coupled to the hydroly- 
sis of adenosine-triphosphate is analyzed in detail showing that it depends 
on the reaction Gibbs energy in a non-linear way. Our results unify ther- 
modynamic and kinetic descriptions, thereby opening new perspectives in 

the study of different transport phenomena in biological systems. 
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1 Introduction 

Mitchell's success in using the proton-motive force to explain energy conversion 
in mitochondria (Mitchell, 1961) led to a considerable effort over the next three 
decades to explore the usefulness of nonequilibrium thermodynamics in biology 
(Caplan and Essig, 1983; Westerhoff and van Dam, 1987; Walz, 1990). Even 
if substantial progress was made, it was eventually concluded that the theory 
had major drawbacks. The flux-force relations that were established, were lin- 
ear, and could only describe the dynamics in a rather narrow range near total 
equilibrium. Since many biological reactions have an activation energy barrier 
and therefore have an intrinsically nonlinear dynamic behavior, such an overall 
thermodynamic analysis could not provide a full description of those processes. 
Existing kinetic descriptions, on the other hand, are not able to account for the 
coupling between transport and chemical reaction. That coupling has so far only 
rigorously been established within the framework of Onsager's theory. These 
difficulties, which already were pointed out by Prigoginc (Prigoginc, 1955) and 
Eyring (Eyring and Eyring, 1965) in the context of activated processes, have 
hampered a wider application of nonequilibrium thermodynamics to biological 
systems, leaving the complete description of active transport processes to ki- 
netic approaches. Kinetic descriptions take the tight coupling of transport and 
chemical reaction as a starting point and can therefore not address a varying 
degree of coupling or slippage. 

Over the last few years, amazing new details have been revealed about the 
molecular mechanism of active transport by some ATPases (Berman, 2001). 
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Experimental evidence has been found that ATP hydrolysis leads to a rotation 
of the central parts of the enzyme, followed by an electric signal in the channel 
that spans the membrane, attributed to the ion that is being transported, see 
e.g. (Peinelt and ApcU, 2004; Burzik et al., 2003; Sun et al, 2004). Active 
transport is accomplished by various types of molecular motors. A theory that 
relates these fluxes with driving forces in the system, however, is still lacking 
(Berman, 2001). 

The purpose of this paper is to show that when irreversible processes are 
analyzed at a mesoscopic level, application of the method of nonequilibrium 
thermodynamics leads to a complete description of the process even in the non- 
linear domain. The mesoscopic level has shorter time and length scales than we 
encounter on the macroscopic level. It is the level of description where fluctu- 
ations matter. Rate laws on this level are stochastic of nature. We shall use 
mesoscopic nonequilibrium thermodynamics (Reguera and Rubi, 2001; Vilar 
and Rubi, 2001) and go beyond the linear overall or 'black box' description of- 
fered by classical nonequilibrium thermodynamics. This gives a new perspective 
to the problem of biological transport phenomena. 

The paper is organized as follows. In section 2, we introduce the main fea- 
tures of active transport of ions and discuss the limitations of a nonequilibrium 
thermodynamic description of the process. In Section 3, we present a kinetic 
description based on thermodynamic grounds which enables us to analyze ac- 
tive transport processes beyond the linear domain. Finally, in Section 4 we 
summarize our main results. 
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2 Active transport and nonequilibrium thermo- 
dynamics 

We shall focus on the problem of active transport of an ion, the meaning of 
nonlinear flux-force relations and the coupling of vectorial and scalar forces in 
this context. As example we take Ca^+ transport in sarcoplasmic reticulum 
by means of the Ca^+-ATPase. We limit ourselves to the conversion between 
chemical energy and transport of the ion, neglecting a possible intervention of 
rotational energy. This does not exclude a torsional mechanism, however (Jain 
et al., 2004). The membrane is highly asymmetric, about 5 nm thick, and 
consists of a phospholipid bilayer with proteins embedded. Figure 1 gives a 
schematic picture of the Ca^^-ATPase embedded in the phospholipid bilayer. 
The reaction takes place at the headgroup of the Ca^+-ATPase. The ion, to 
be actively transported along the protein channel, binds to a separate site and 
triggers the chemical reaction, the ATP hydrolysis. Each channel transports 
only a few molecules or ions per second, and there are not many channels of 
the same type in one membrane. Experiments arc done in a solution with many 
vesicles or organelles, so the observation refers to the average over the channels 
(Burzik et al., 2003). Thus, we have a chemical reaction on a time scale similar 
to that of diffusion (Peinelt and Apell, 2004; Burzik et al., 2003). Can we 
legitimately consider this dynamic situation with thermodynamic tools? We 
shall see that the answer is yes. 

The spontaneous chemical reaction taking place at the membrane surface in 
Fig.l is, schematically: 

ATP (s) ^ ADP (s) + P (s) 
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Figure 1: The Ca^+-ATPase that transports Ca^+ by hydrolysis of adenosine 

triphosphate (ATP) to adenosine diphosphate (ADP) and phosphate (P). The 
reaction takes place at the enzyme headgroup, and triggers the active transport 
through the protein stalk and channel. A surface area element is indicated by 
the dotted line. 

where s stands for surface. The reaction Gibbs energy which arises is 



where jX^ is the chemical potential at the surface of component j . Calcium ions 
are transported across the membrane from side i to side o. The transport takes 
place against the chemical potential gradient as a consequence of its coupling 
to the above reaction: 



There is probably also a countertransport of protons to keep the system's local 
electroneutrality (Peinelt and Apell, 2004). 

Non-equilibrium thermodynamics provides a description of the irreversible 
processes occurring at the surface (Albano and Bedeaux, 1987). It has been 
shown that the surface behaves as an autonomous thermodynamic system in 



AG — jlp + //.ADP ~ A*ATP 



(1) 



Ca2+ (i) ^ Ca?+ (s) ^ Ca?+ (o) 
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local thermodynamic equilibrium, even in the presence of strong driving forces 
(R0sjorde et al., 2000). This fact enables one to formulate the Gibbs equation 
for the surface element indicated in Fig.l (Albano and Bedeaux, 1987): 

Tds^ — du^ — A«ATP'^c5^TP ~ Madp^'^adp ~ fJ-pdcp ~ A^Ca^C^a (2) 

Here c| is the surface concentration of component j in mol/mg, referring to the 
total unknown area of transfer per mg of material, and u*^ is the internal energy 
density of the surface. If the membrane potential contributes to transport of 
the calcium ions, one has to use electrochemical potentials. The conservation 
laws for the reacting species in Eq.lQ) are 

^^±. + j;. (3) 

where Jj are fluxes of species, r is the reaction rate (all in mol/s mg), and t is 
the time. The minus sign refers to ATP and the plus to the other species. The 
concentration of Ca^+ in mol/mg on the surface satisfies 

dcce. _ ji 



dt 



- Jh. ~ Jc. (4) 



where are the inwards and outwards ion fluxes in mol/s mg. We omit the 
charge number in subscripts. 

For the time being, we consider an isothermal surface. The entropy produc- 
tion (in W/mg K) then follows by using Eqs.lO and l@J in combination with 
Gibbs equation, taking the temperature constant and /i| = /ij : 

a^ = -rlAG- JSai(/.°ca-Mk) (5) 

The countertransport of protons to maintain the system's electroneutrality (Peinelt 
and Apell, 2004) can be accounted for by modifying the driven force for Ca^+, 
see Eq. (20) below. From this one normally has inferred linear laws between 



fluxes and driving forces (Caplan and Essig, 1983), with an expression for the 
flux of Ca2+ 

JS, = -:^AG-:^(M°ca-A^ba) (6) 

describing hnear active transport for small values of AG and constant Onsager 
coefficients Ldr, Lrr- Equation (jSJ is normally taken for a 'black box' description 
of active transport fromside i to o (Caplan and Essig, 1983). When the second 
term is negative, the coupling of Jq^ to the reaction (via Ldr) can compensate 
that negative contribution thus making transport against the concentration gra- 
dient possible. This is the coupling phenomenon that is referred to as active 
transport in the literature. No explanation has yet been given for the nature of 
the coupling (Berman, 2001), leaving the problem unsolved as how to charac- 
terize the total process far from equilibrium. 



3 Kinetic description from thermodynamic grounds 

The lack of equilibrium for the chemical reaction on a time scale similar to that 
for ion transport is the reason for energy transduction between reaction and 
transport. This energy transduction process shall now be described for shorter 
time scales at which the system evolves from the initial to the final states. The 
local states of the system correspond to successive molecular configurations and 
are parametrized by an internal reaction coordinate 7 running from 0, at which 
there are reactants, to 1, at which there are products. The coordinate 7 is 
dimensionless. One may then extend the assumption of local thermodynamic 
equilibrium to the mesoscopic level and formulate the corresponding Gibbs equa- 
tion in 7-space: 

Tds^(7) = -G(7)dc.(7) - M^a(7)rf4a(7) (7) 
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where 0^(7) is the concentration of enzymes with a reaction mixture charac- 
terized by coordinate 7, and 0^3^(7) is the average concentration of Ca^+ over 
the state given by 7, both concentrations in mol/mg. Before the reaction, the 
Gibbs energy becomes G(0)y(iATP ; after the reaction it is 0(1) = /xp + /xadp- 
The relation between AG and 6(7) is shown to the left in Fig. 2. 




1 y, 1 

Y Y 

Figure 2: Left, the expansion of 6(7) in reaction coordinate space. Right, the 
activation energy $(7) with a peak at 74^- 



Mass conservation can be expressed at any position in 7-space through the 
conservation laws 

^_|.,(,) ^ = Jk(7)-JSj7), (8) 

Here r{'y) is the local reaction flux along the reaction coordinate. The fluxes 
Jcail) and Jca.il) contributions to the diffusion fluxes on sides i and o respec- 
tively, averaged over enzymes in state 7. The integrals of the fluxes over 7 give 
the corresponding global fluxes. This choice for the conservation equation for 
Ca^+ may capture the essence of the delivery of Ca^+ to the other side in a 
pump that allows slippage. 

The local entropy production ct*(7), defined through the entropy change per 
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unit of time in the isothermal surface is then: 

— ^^=Cr(7) + - ^^^^ ^ —Jcail) - —Jcah)^ (9) 

then foUows from Eqs.Q and ^ 

^^(7) = -K7)^^^ - ^ JSa(7)(M°Ca - ML(7)) - ^ Jca(7)(ML(7) - Mba) 

(10) 

As has been assumed aheady, there is equihbrium for the Ca^+-ion between the 
surface and the bulk solution /ica(7) ^ I^Ca- "^^^ local entropy production then 
simplifies to 

^^^(7) = -K7)^^ - ^JcMi^^h. - i^U (11) 

Due to the activation barrier, the reaction rate almost immediately reaches a 
quasi-stationary state in which it only depends on time. Using this property, 
the entropy production (4) is recovered upon integration of Eq. with a 

varying calcium flux. 

We shall assume that Ca^"*" is not transported across an activation energy 
barrier. This means that the ion can escape from any surface state given by 7 
with /ica(7) ~ Mca to side o with ij-q^^, an explanation for why the flux depends 
on 7. The local flux- force relations that follow from Ea. (|ll|l . are: 

rh) - -^^-^(M°Ca-.ba) 
^Sa(7) ^ -^^-^(.°Ca-.ba) (12) 

which are valid along the path going from the initial to the final state and have 
been formulated under the equilibriation assumption discussed previously. At 
each point of the path, the Onsager relation Irdil) = Idril) is fulfilled and the 
coefficient do not depend on the driving forces. The last relation will give a 
non-linear equation for the flux of Ca^+ across the membrane. 



The chemical reaction is an activated process. Along the internal coordinate, 
we can write the Gibbs energy as a combination of the ideal contribution and 
an activation energy per mole, <&(7); 

G(7) = G" + i?r In c(7) + $(7) (13) 

The boundary values of £'(7) were defined above. The standard Gibbs energy is 
chosen to be G° = Matp- I^^r the potential profile, correct boundary conditions 
are obtained with $(0) = and $(1) = /^adp + Mp- The activation energy 
barrier with the transition state, "ftr, is illustrated in Fig. 2. 

We do not know the coefficients in Eas. (|12|) . But, under the condition of 
constant stoichiometry, or tight coupling, we can make certain assumptions that 
allow us to integrate the expression for the calcium flux. A constant stoichiome- 
try means that the number of Ca^+ transported for every reacted ATP is fixed. 
When the stoichiometry is fixed, Jq^ and r are dependent, and the matrix of 
coefficients has a zero determinant. The number of Ca^+ transported for every 
reacted ATP can also be defined in 7-space. We obtain: 

n(7) = ^H4 = rTT (1^) 

r lrr(J) Irdyi) 

In order to show how the linear, local Eas. H12|l transform into nonlinear global 
laws, we choose as thermodynamic force the gradient of the local fugacity z{j), 

z(7) = exp[(G(7) - G")/i?T] = 0(7) exp[$(7)/i?T] (15) 

The fugacity divided by the fugacity at equilibrium, can be understood as the 
probability distribution function in 7-space divided by the equilibrium distri- 
bution, Cr(7)/cr.eg. The usc of probabilities is characteristic to a mesoscopic 
description. The local diffusion fiux (|12|l can then be rewritten in terms of the 
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new thermodynamic force 

TO / ^ raD^ dz{-i) Iddil) r o i 1 cifiA 

-^c-^^'- exp[ci>(7)/i?r] 97 ^^L'^ca McaJ (ibj 

where we have used Eas. (|14ll and (|15|) and defined D' — lrr{'j)R/c{'j) in accor- 
dance with Kramers [?]. We find the macroscopic flux in terms of the driving 
forces by integrating over 7. Using the quasi-stationary condition for the re- 
action flux, and constant n, it follows that Jg^ is constant. Using furthermore 
that D' is in good approximation constant, we obtain: 

^Ca = -nDcATP ( exp ^ - 1 ) - Lca [/i&a ~ MCa] (17) 



where 



D = D\[ exp$/i?r)-i and ica = / Idd/Tdj (18) 

JQ Jo 

We have also used exp(G(0) - G°)/RT = catp- 

This overall diffusion flux H21|l describes active transport of Ca?^ for arbi- 
trary values of AG, and accounts thus for the transduction process also in the 
nonlinear regime, which could not previously be described through nonequilib- 
rium thermodynamics. Nonlinearities inherent to the kinetics of the chemical 
reaction emerge in the broad description if one proceeds with a thermodynamic 
analysis at the mesoscopic level and subsequently coarsens the description by 
retaining only the initial and final states. At small values of the driving force. 
Eg. 1)21(1 reduces to expression (5) valid only in the linear regime. It is interesting 
to note that the prefactor of the nonlinear term contains an Arrhenius factor 
characteristic for activated processes. The first term in Ea. (|21|l is also present 
in the equation for active transport, derived by the King-Altmann-Hill method 
(Westerhoff and van Dam, 1987). This equation does not have a second term. 
The derivation and the understanding of Eg. ((211) compared to equations that 
have been used earlier to describe Ca^+ transport, however, are very different. 
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Equation 12 1|) can be used to plot experimental data and find the product 
nD and £ca- The value of n obtained by imposing the static head condition 
"^Ca = to Eq. (17), may be different from the one given through the relation 
n = J^a/r depending on experimental conditions. The chemical potentials for 
Ca^+ for state i or o: 

M^°a = /^Ca+^rlnC^: (19) 

where the ideal mixture approximation and the same standard state has been 
used for both sides. The approximate difference in chemical potentials for Ca^^ 
is: 

M&a - Mba - i?T In % ^ £Ca_£Ca (2O) 
^Ca ^Ca 

This driving force can be modified to also account for the countcrtransport of 
protons. The effective driving force for exchange of Ca^+ and H+ becomes 

Afi = RT In (21) 

Only a zero difference in pH between the sides has no influence on this force. 

In the derivation of the flux equations H12|) . we have made the assumption 
that the surface is isothermal. Otherwise the model is general in the way that 
it allows slippage. We have further assumed that Ca^+ transport across the 
membrane is not an activated process, but rather that its contribution to the 
total flux is simply linear in its corresponding driving force, with a transport 
coefficient independent from the activation energy. Therefore, the only nonlin- 
earity comes from the kinetics of the chemical reaction. This assumption may 
be valid for the concentrations of calcium that are physiologically relevant. 

In order to integrate the flux on the mesoscopic level, we introduced one 
more assumption, namely that of a tight pump. This assumption is thus a 
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requirement for Ea. H21|) . The assumption may hold for an isothermal surface, 
but there are many conditions for which this is not true (Berman, 2001). 

The theory we have presented could be useful in describing other types of 
active transport in biological systems, as the ones taking place in molecular 
motors for which translational (Bedeaux et al. , 2004) or rotational motion may 
be induced at the expense of the energy provided by the chemical reaction (Sun 
et al., 2004). It can also be used for a more detailed look at pump slippage, be- 
lieved to occur for high driving forces, and indicated by a variable stoichiometry 
(Zoratti et al., 1986; Walz, 1990). One advantage of the present formulation, not 
contained in kinetic approaches, is that the energy dissipated as heat (Berman; 
2001) can be directly calculated from Ta^ . A large value of Ta^ suggests that 
the surface may be nonisothermal, and Ta^ has a meaning also for nonlinear 
flux- force relations. 

4 Conclusions 

The main aim of this paper has been to show that active transport processes can 
be described by means of mesoscopic nonequilibrium thermodynamics (Reguera 
and Rubi, 2001; Vilar and Rubi, 2001). We have been able to describe the 
coupling between the chemical reaction and transport, which is mandatory in 
a complete thermodynamic description of this type of system, and which was 
previously not obtained outside the range of linear flux-forcc-relations (Berman, 
2001). 

We have shown that when local equilibrium is formulated at shorter time 
scales, the overall diffusion flux is nonlinear in the macroscopic driving force. A 
new way to plot experimental data has thus been given, applicable for large 
values of AG. The formulation we propose is of a mesoscopic nature and 
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may incorporate the presence of fluctuations in the dynamics consistent with a 
Fokker-Planck description (Vilar and Rubi, 2001). We have somewhat arbitrar- 
ily taken active transport of Ca^+ as an example, but the scheme proposed could 
be transferable to many kinds of transport processes that obtain their energy 
from chemical reactions such as molecular motors. This shows that nonequilib- 
rium thermodynamic methods are not restricted to the linear response domain 
(Jiilicher et al., 1997) but can be used in a broader context to perform a de- 
scription of the dyamics also in the nonlinear domain. We hope in this way to 
provide a general scenario in which nonlinear transport processes occurring in 
biological system can be dealt with. 
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